Coupled mode theory for photonic band-gap inhibition of spatial instabilities 
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We study the inhibition of pattern formation in nonlinear optical systems using intracavity pho- 
tonic crystals. We consider mean field models for single and doubly degenerate optical parametric 
oscillators. Analytical expressions for the new (higher) modulational thresholds and the size of the 
"band-gap" as function of the system and photonic crystal parameters are obtained via a coupled- 
mode theory. Then, by means of a nonlinear analysis, we derive amplitude equations for the unstable 
modes and find the stationary solutions above threshold. The form of the unstable mode is different 
in the lower and upper part of the bandgap. In each part there is bistability between two spatially 
shifted patterns. In large systems stable wall defects between the two solutions are formed and we 
provide analytical expressions for their shape. The analytical results are favorably compared with 
results obtained from the full system equations. Inhibition of pattern formation can be used to 
spatially control signal generation in the transverse plane. 



I. INTRODUCTION 



Photonic crystals (PC) have shown to be able to control light in ways that were not possible with conventional optics 
P, Q . Their remarkable properties stem from the unusual dispersion relation as a result of the periodic modulation 
of their dielectric properties. Most of the interest on PC is related to propagation problems. Here the existence of 
photonic band-gaps, i.e. a range of frequencies for which light can not propagate in the medium, allows for a new 
way of guiding and localize light once defects are introduced in the PC Q, S- Photonic crystals in combination with 
— h ■ nonlinear effects have been also considered for all-optical switching devices Q- 

Transverse effects in periodic media have been mainly studied in propagation in planar waveguides with periodic 
modulation of the refractive index in the transverse direction and arrays of couple waveguides. In the first case the 
attention has been focused on the so called spatial gap (or Bragg) solitons 0, S El > which are intense peaks with 
frequencies inside the (linear) band-gap. This is possible thanks to a shift of the photonic band-gap boundaries due to 
nonlinear effects: the intense core of a gap soliton can propagate freely in the periodic media, while in its less intense 
tails, where non-linearity can be neglected, light is reflected back to the center due to Bragg reflection, sustaining the 
localized structure. Bragg solitons are usually studied within the coupled-mode theory where only the slowly varying 
envelope of two counter-propagating beams are considered. Arrays of waveguides are usually studied within the tight- 
ly ■ binding approximation, where the system is described by a set of coupled ordinary differential equations describing 
the dynamics of the guided-mode amplitude in each site. The evanescent coupling between adjacent waveguide give 
rise to the so called discrete diffraction. If nonlinearity is also present discrete solitons can be formed 0, 0] • 

In this paper we consider a different case, a photonic crystal inside a nonlinear optical cavity, i.e. a system with 
driving and dissipation. Nonlinear optical cavities typically undergo spatial instabilities leading to the formation of 
spatial structures similar to what observed in many different fields across science In particular, spatial structures 

.{pH \ in nonlinear optical cavities have important potential applications in photonics such as memories, multiplexing, optical 
processing and imaging Control of spatial structures is then an important issue for the implementation of such 
devices. In this context, different mechanisms of control have been proposed 0,0,0]- In particular, Neubecker and 
Zimmerman reported on experimental pattern formation in presence of an external modulated forcing and observed 
lockings between the forcing and natural wavelengths 0] . More recently the use of the properties of photonic crystals 
for the control of optical spatial structures has been considered [T^ IT'S. ITU . I n 0] we showed how the photonic band- 
gap of a photonic crystal could be used to inhibit a pattern forming instability in a self-focusing Kerr cavity. Pattern 
formation in an optically-injected, photonic-cr ysta l vertical-cavity surface-emitting laser, electrically biased below 
threshold has been observed experimentally in [14(. Finally, in [15|, a defect in a photonic crystal has been used to 
achieve stable emission in broad area lasers. 

Here we study the inhibition of pattern formation in a mean field model for a Degenerate Optical Parametric 
Oscillator with a photonic crystal. We first show that the inhibition mechanism introduced in ^| occurs independently 
of the type of nonlinearity here being quadratic. Then we introduce an analytical treatment by means of a couple- 
mode theory that fully explains the phenomenon. A linear stability analysis allows us to determine analytically 
the new thresholds shifted in parameter space by pattern inhibition and the form of the unstable modes. We show 
that the band-gap region is divided into two different regions: one in which pattern formation takes place due to 
an energy concentration at the maxima of the photonic crystal modulation, and a different one where the energy 
concentrates at the minima. We determine amplitude equations for the unstable modes and find the stationary 
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FIG. 1: Scheme of a doubly resonant degenerate optical parametric oscillator with a photonic crystal. It consist of a ring 
cavity filled with a quadratic (x^) medium and a photonic crystal slab (PC). E is the plane wave input field at frequency lo 
partially transmitted into the cavity. The other mirrors are assumed to be perfectly reflecting. B and A are the pump and 
signal fields at frequencies uj and u/2 respectively. In the single resonant case only the signal field A resonates with the cavity. 



solutions above threshold analytically. We show that this scenario is organized by a co-dimension two bifurcation 
point where both modes become simultaneously unstable. We also show that in the band-gap, where the pattern 
arises with a wavenumber half the one of the PC, there is bistability between two different "phase locked" solutions 
differing by a shift of a PC wavelength in the near field. This bistability stems from the breaking of the translational 
symmetry by the PC. In large systems this leads to the formation of domains of different "phases" connected by 
topological walls (or antiphase boundaries) Q . This is a different kind of (spatial) bistability where the two different 
"phases" do not differ on the phase of the electric field but on the location of the spatial oscillation in the near field. 



II. MODELS 

A phase-matched Double Resonant Degenerate Optical Parametric Oscillator (DRDOPO), where both pump and 
signal fields resonate in the cavity, can be described in the mean field approximation by |16| : 

d t B = T[-B + E-A 2 ]+i^d xx B 

d t A = —A — iA s A + BA* + id xx A, (1) 

where B and A are the pump and signal slowly varying amplitudes, T is the ratio between the pump and signal 
cavity decay rates and E is the amplitude of the external pump field (our control parameter). A s = A + a fix) is the 
total signal detuning, where A is the average detuning and f(x) = e lkpaX + e~ lk *' aX describes the spatially dependent 
contribution of the photonic crystal (Fig. ^) . 

For the singly resonant DOPO (SRDOPO), where there is no cavity for the pump field, the mean field equation for 
the resonating signal is 0, Il8j 

d t A = -A-iA s A + EA* - \A\ 2 A + id xx A. (2) 



III. LINEAR STABILITY ANALYSIS WITH PERIODIC MEDIA: COUPLE-MODE THEORY 

The linearization around the steady state homogeneous solution B° = E and ^4° = of Eqs. (JTJ and J5J leads, in 
both cases, to the same equation for the perturbations of the signal A: 

d t A = -(1 + iA)A - iaf(x)A + EA* + d xx A. (3) 

Without the photonic crystal (a — 0), the homogeneous solution is stable for E < 1. Above threshold (E > 1) a 
stripe pattern arises with a wavenumber k c = \J — A. Here we restrict ourselves to negative values of the detuning. 
For positive detuning the system display bistability between homogeneous solutions 18]. For a ^ 0, and assuming 
that the amplitude of the modulation is weak enough, we can write the signal perturbations as a superposition of two 
waves with opposite transverse wavenumbers 



A = A+{x, t)e lk "" x/2 + A_(x, t)e~ ik ^ x/2 , 



(4) 
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where A± (x, y) are slow functions of x. This is basically equivalent to the so called coupled-mode theory for propagation 
of pulses in periodic media [3J. 

By writing A± = a^{t)e ikx + ag(t)e- lkx , one obtains a set of coupled linear ordinary differential equations for 
the amplitudes a± (a±) of the Fourier components ±fc pc /2 + k (±k pc /2 — k) of the perturbations. By defining 
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this set of equations can be written as 

F = £F, 



where C is given by 
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In this way one reduces the stability analysis of the steady state of the initial partial differential equations with 
periodic coefficients and (0) to diagonalize the 4x4 complex matrix C. In the case k = one has to consider that 



The eigenvalue of C with largest real part is 
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and by setting A = we obtain four marginal stability curves 

A^k.E) = -dt(k) + d 2 {k) + d 3 {E) 
A 2 (k,E) = -di(fc) + d 2 (k) - d 3 (E) 
A 3 (k,E) = -d 1 (k)-d 2 (k)+d 3 (E) 
Ai(k,E) = -di(fc) - d 2 (k) - d 3 (E). 



(8) 



where d x {k) = (k pc /2y + k 2 , d 2 (k) = yj Ak 2 (k pc /2) 2 + a 2 and d 3 (E) = y<E 2 - 1. 

Fig-Elshows the marginal stability curves for different values of the pump E. Note that, due to definition 0J, k = 
indicates a perturbation with a wavenumber at the limit of the first Brillouin zone. Therefore it is convenient to plot 
A as a function of k' — (k pc /2 — k)/k pc . Dotted lines are the results from the coupled-mode theory (JHJ), while solid 
lines have been obtained from a numerical stability analysis of the full model, i.e. solving the eigenvalue problem 
associated to the linear differential operator with periodic coefficients in the rhs of Q 13]. The coupled- mode theory 
provides a very good analytical approximation for thresholds and unstable wavenumbers allowing us to predict the 
existence and size of a band-gap in the modulation instability analytically. In the following we analyze the results of 
the coupled-mode theory in more detail. 

From © one can evaluate the instability threshold of the fundamental solution as a function of the detuning A. For 
E < 1.0 the steady state homogeneous solution is stable for all values of the detuning. At E = 1 (Fig.^) d 3 (E) = 
and the four marginal stability lines given by Eqs. JHJ become only two, signaling the instability threshold for all 
values of the detuning outside the band-gap (— dx(k pc /2) — d2(k pc /2) + a 2 < A < —(k pc /2) 2 — a < —(k pc /2) 2 + a < 
A < —d\{k pc /2) + d 2 {k pc /2) + a 2 ). In the bandgap, —(k pc /2) 2 — a < A < —(k pc /2) 2 + a, the instability is, however, 
inhibited. The threshold E = 1 for values of A outside the bandgap is slightly underestimated. In the full model, 
for this value of the pump, the system is still stable (no solid line in Fig[2J0- This is due to the fact that the spatial 
modulation couples the fundamental wavenumbers with their harmonics, which are dumped, introducing an additional 
source of stability. Harmonics are not taken into account in the couple-mode theory and then the threshold is slightly 
lower than for the full models Q and J2J). 

Increasing further the value of the pump the band-gap narrows (Figs |2d,c). At E — \Jl + a 2 the homogeneous 
solution become eventually unstable for any value of the detuning (Fig.|2Jl). The threshold value of E as a function 
of the photonic crystal parameters for any value of A is given by 
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FIG. 2: Unstable wavenumbers of a DOPO for a)E = 1.0, b)E = 1.02, c)E = 1.05, d)E = Vl + a 2 = 1.118034, and e)E = 1.2 
in presence of a periodically modulated media (k pc = 2, a — 0.5). For each value of A wavenumbers between the two lines are 
unstable. 
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FIG. 3: Pump threshold as function of the detuning from the theoretical result @ obtained by means of the coupled-mode 
theory (solid black and grey lines). The dashed line is the result obtained from a numerical stability analysis of the full model. 
The dotted line is the threshold for a DOPO without photonic crystal. 

and the critical wavenumber by: 



(fc pc /2) 2 -A- ^-4(^)2A + a2 -{^f + a <A< -di(^f) + d 2 (^) 

-(%^) 2 <A< _(^)2 + a 

-{^) 2 -a <A< -{h^f 

(i£)2_ A _^_ 4 (i£)2 A + a 2 _ dl (^)2_ d2 (^ £)+a 2 <A < _(^££)2„ a 



(10) 



Fig. El shows the theoretical prediction for the instability threshold of the homogeneous solution (solid lines). 
The black and grey lines correspond to two different eigenmodes. The first is a mode with energy concentrated in the 
minima of the photonic crystal, while the second in the maxima (see next section for more details). The dashed line 
is the threshold for the full system computed as explained above. 

The coupled-mode theory accurately capture the features of the inhibition of the modulation instability by the 
photonic crystal. The fact that the gap of values of the detuning for which pattern formation is inhibited decreases 
with the intensity of the pump is because nonlinearity always overcomes the linear inhibition by the photonic crystal 
for suitable high input intensities. This is similar to the phenomenon behind the formation of gap solitons. But 
in that case the nonlinearity only shifts the position of the band-gap, while in our case it narrows it and makes it 
eventually to disappear. 

From Eq. © we can calculate the size g(E, a) of the band-gap as function of the system and photonic crystal 
parameters: 



g(E,a) = 2(a~VE rZ l)- (11) 

The size of the band-gap is proportional to the amplitude of the modulation a and becomes smaller by increasing the 
pump E (Fig. [21, eventually disappearing for E = VT + c? (Fig. EJi). For E > 1 + a 2 the homogeneous solution is 
unstable for any value of A (Fig. [2k) . Fig. 0] shows how the width of the band-gap (white region) changes as function 
of the amplitude of the modulation a for a fixed value of the pump E. The presence of the modulated medium 
opens an entirely new stable region. Fig. 0] should be compared with Fig. 6 of [l3| showing that pattern inhibition is 
independent of the form of the nonlinearity. 

Since for negative signal detunings, in the absence of PC, down conversion takes place at a finite wavenumber, the 
inhibition mechanism can be used to spatially control the generation of signal in a way analogous to that explained 
in [l3|. If we set our system in a parameter region where pattern formation is inhibited, the inclusion of a defect in 
the PC will lead to a spot of signal generation as shown in Fig. 
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FIG. 4: Stability diagram of the steady state homogeneous solution for E = 1.1 (solid line). The system is unstable in the 
shadowed region. The dashed (dotted) line shows the boundary of the stable region for E = 1.0 (E = 1.2). 
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FIG. 5: Localized spot of signal generation due to a defect in the photonic crystal. Here E = 1.08, A = —l,kp C = 2 and 
a = 0.5. 



IV. WEAKLY NONLINEAR ANALYSIS: AMPLITUDE EQUATIONS FOR PATTERN FORMATION IN 

PERIODIC MEDIA 

In this section we study by means of a multiple scales analysis the solutions that appear above threshold for values 
of the detuning inside the band-gap. Above threshold nonlinear terms have to be considered since they saturate the 
linear growth induced by the instability. In the following we consider the SRDOPO only. By including the nonlinear 
terms from J5J in © we obtain: 

F = £F-W(F), (12) 

where W = [(a^a^* + 2aZ*al)a+, (2a^al* + a^* a+)aZ* (a'taZ* + 2a+*at)a^, (2ata+* + aC*a^)a7*] T is a nonlinear 
function of F. 

In the bandgap, the critical wave number is always k c = independently of the detuning. We recall that in this 
case one has to consider a+ = a~. In the following we will discuss the lower (— (fc pc /2) 2 — a < A < —(k pc /2) 2 ) and 
upper (—(kp C /2) 2 < A < — (k pc /2) 2 + a) halves of the band-gap separately. 

A. Lower-half part 

Assuming the following scaling |l6j : 





F = eFi+e 3 F 3 

E = E t h + e 2 i?2 



(13) 
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T = e 2 t 
X = e 2 x. 

Substituting (|13|l in l|12fl . at order e one obtains Fi = FiVi, where vi is the critical eigenmode of C 
(£vi = and |vi| = 1) and F\ is its real amplitude. In the near field, vi has the form: [(1 + E t h — 
i^/Ej h - l)/y/2Et h (l + E th )] cos(fc pc x/2). 

The solvability condition at order e 3 yields to the following equation for the real amplitude A\ — eFi of the unstable 
mode 

d t A x = y/E* h - 1 d 2 x A x + E th (E - E th )A 1 - ^A* . (14) 

The amplitude equation (|14fl is equivalent to the one obtained from a secondary instability at twice the spatial period 
of a cellular pattern |19| . except for the fact that in our case the translational invariance has been broken by the 
presence of the photonic crystal. The homogeneous steady state solution of (|T4"|) is A\ = ±2^ E t h(E — E t h)/3. The 
solution of J2J is then: 

a{x) = ± v!vtt5 + Eth - - i ) cos{k -f^- ( i5 ) 

The bifurcation diagram and spatial form of this solution is shown in Fig[B] The plus and minus solutions l|15f) are 
created in a pitchfork bifurcation. The analytical solution (|15f) is in very good agreement with the stationary solution 
of the full model computed numerically. Note that, despite being completely equivalent, the plus and minus solutions 
of (|15fl are not the same. One corresponds to the other shifted by a wavelength of the photonic crystal modulation. 
In a system with translational invariance the position of a solution that brakes such symmetry is undetermined. The 
photonic crystal periodicity selects just two of a continuum of possible solutions. This two " frequency locked" solutions 
differ now by a shift of a photonic crystal wavelength in the transverse position. This situation is the spatial analogue 
of an oscillatory system forced at twice its natural frequency |20j . In a large system and starting from arbitrary initial 
conditions, some sections of the system will attain the first solution while others will move to the second one, leading 
to the formation of domain walls. This is illustrated in Fig. [7] The solid line is the result of a simulation of the full 
model @ starting from a random initial condition. The dashed line is the steady state solution of Eq. (|14fl connecting 
the two (plus and minus) homogeneous solutions |8| times the critical eigenmode vi in the near field: 

Mx) = vfVfrlir (* + Eth i v / ^ TT ) tanh 

The dotted-dashed line in Fig. is the envelope of (|16|l . The analytical result (|lf>|) is in very good agreement with 
the domain wall obtained from the numerical simulations of the full model J3J. 
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B. Upper- half part 



In the upper-half part of the photonic band-gap the critical mode vi has the form: {[y/Ef h — 1 + i[Eth ~ 
\)\j \j2E t h{Eth — 1)} sin(fc pc x/2). As in the previous case, we obtain a similar amplitude equation for the real ampli- 
tude A\ of the unstable mode 



dtA x = \ E 2 th - 1 dlAx + y/E - EthAx 



3 ~ 3 



(17) 



In this case the pattern solutions is 



E — E t h 
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E th 
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(18) 



The bifurcation diagram and spatial form of this solution is shown in Fig EI The analytical solution (|18|l is in very 
good agreement with the stationary solution of the full model computed numerically. As in the previous case, in large 
systems, fronts between the plus and minus solutions are formed. The shape of the front is given by: 



A x {x) = 2 
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3V E t h(E th -l) 
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FIG. 6: Left: Bifurcation diagram for the pattern arising above threshold in the lower-half of the band-gap. Here A = —1.2. 
For this value of the detuning Eth = 1.04403. Right: Real (solid line) and imaginary (dashed lined) of the pattern solution of 
the full model for E = 1.2. The dot-dashed and dotted lines correspond to the real and imaginary parts of the pattern from 
the coupled-mode theory. The grey solid line illustrates the modulation of the photonic-crystal for comparison. 
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FIG. 7: Two domains corresponding to the plus and minus sign solutions 1150 separated by domain walls. The final state is 
the result of a numerical simulation starting from an arbitrary initial condition. Left: Real part of the field. Right: Close up to 
the intensity around a domain wall. The dotted line shows the modulation of the photonic crystal. Here A = —1.1, E = 1.15, 
a = 0.5 and k pc — 2.0. 



Note that while the cosine solution in the lower-half part of the band-gap has the maxima of the intensity at the 
maxima of the photonic crystal modulation, the sine solution in the upper-part has them at the minima. This different 
distribution of energy in the photonic crystal is at the basis of the creation of the band-gap. An interesting point is 
the middle of the band-gap (A = — (k pc /2) 2 ). At this particular value of the detuning, for E = E th = vTTe? both 
the sine and cosine modes become simultaneously unstable. This is a co-dimension two point where Eqs. I|14fl an d 
(|17ll become coupled. The unfolding of such a critical point is however beyond the scope of this paper and it is left 
for future investigation. 



V. CONCLUSIONS 



In this paper we have studied pattern formation in nonlinear optical cavities in presence of a photonic crystal, i.e. 
a spatial modulation of the refractive index. The linear phenomenon of the band-gap inhibits pattern formation for a 
certain range of cavity detunings (bandgap). For high enough intensities nonlincarity finally overcomes the inhibition 
by the photonic crystal and a pattern arises. By means of a couple- mode theory approach we have obtained analytical 
expressions for the new (shifted) threshold and the form of the unstable modes. The bandgap is naturally divided 
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FIG. 8: Left: Bifurcation diagram for the pattern arising above threshold in the lower-half of the band-gap. Here A = —0.9. 
For this value of the detuning E t h = 1.07703. Right: Real (solid line) and imaginary (dashed lined) of the pattern solution of 
the full model for E = 1.2. The dot-dashed and dotted lines correspond to the real and imaginary parts of the pattern from 
the coupled-mode theory. The grey solid line illustrates the modulation of the photonic-crystal for comparison. 



in two halves, the lower-half one in which the unstable mode has a cosine shape, i.e. its intensity maximums are in 
correspondence with the maximums of the photonic crystal modulation, and the upper-half where the unstable mode 
is a sine; the intensity maximums are in correspondence with the minima of the photonic crystal modulation. By 
means of a multiple scale analysis we also found the pattern solution above threshold. In each part of the bandgap 
there is bistability between the plus and minus cosine or sine solutions. This bistability stems from the breaking of the 
translational symmetry of the photonic crystal. In large systems domain wall between this two solutions are typically 
form. The shape of the defect wall is given by an hyperbolic tangent. While the particular form of the coefficients 
are model dependent, the shape of the unstable modes and the splitting of the bandgap in two different regions are 
generic. We have checked that the same phenomenon is present in a completely different model, namely the Kerr 
cavity model studied in [l3|. Finally we also have shown that photonic crystal can be useful to engineer particular 
spatial signal outputs in frequency down conversion. An interesting extension of this work is to consider the case with 
two transverse dimension. In this case, some of the features studied here, such us the two different modes in each part 
of the bandgap, will remain the same, while new features like the coupling of the geometry of the photonic crystal 
and that of the spontaneous pattern will come into play. 
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